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A generic second-order scalar-tensor theory contains a nonlinear derivative self-interaction of 
the scalar degree of freedom (p a la Galileon models, which allows for the Vainshtein screening 
mechanism. We investigate this effect on subhorizon scales in a cosmological background, based on 
the most general second-order scalar-tensor theory. Our analysis takes into account all the relevant 
nonlinear terms and the effect of metric perturbations consistently. We derive an explicit form of 
Newton's constant, which in general is time-dependent and hence is constrained from observations, 
as suggested earlier. It is argued that in the most general case the inverse-square law cannot be 
reproduced on the smallest scales. Some applications of our results are also presented. 
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I. INTRODUCTION 

Since the discovery of the accelerated expansion of the 
Universe numerous attempts have been proposed to 
explain the origin of this biggest mystery in modern cos- 
mology. A vast class of models for cosmic acceleration 
invokes a scalar degree of freedom, (f>, which may cou- 
ple minimally or nonminimally to gravity and ordinary 
matter. In the case of nonminimal coupling, the mod- 
els are commonly called modified gravity (or "dark grav- 
ity") rather than dark energy, as (jj participates in the 
long-range gravitational interactions and thereby accel- 
erates the cosmic expansion. Modified gravity models 
must be designed with care, because otherwise the effect 
of modification could persist down to small scales, which 
could easily be inconsistent with stringent tests in the 
solar system and laboratories. For this reason, screening 
mechanisms for scalar-mediated force are crucial. 

There are mainly two approaches for screening the 
scalar degree of freedom in modified gravity models. The 
first one is the Chameleon mechanism [2|, by which the 
scalar acquires large mass in a high density environment. 
This is employed in viable f{R) gravity 0], which is 
equivalent to a scalar-tensor theory with an appropri- 
ate potential. The second one is the Vainshtein mech- 
anism In this case ^'s kinetic term becomes effec- 
tively large in the vicinity of matter due to some nonlin- 
ear derivative interaction, suppressing the effect of non- 
minimal coupling. The Vainshtein screening is typical 
in Galilcon-likc models Q and nonlinear massive grav- 
ity (e.g., @). In the Vainshtein case, nonlinearities play 
an important role in possible recovery of usual gravity 
on small scales even in a weak gravity regime. To test 
models of modified gravity against experiments and cos- 



mological observations, we therefore need to clarify the 
behavior of gravity around and below the scale at which 
the relevant nonlinearities set in. In this paper, we ex- 
plore the consequences of the latter mechanism in detail, 
taking into account the nonlinear effect. 

We study gravity sourced by a density perturbation 
of nonrelativistic matter, on subhorizon scales in a cos- 
mological background, using the (quasi)static approx- 
imation. To provide generic results, we work in the 
most general scalar-tensor theory with second-order field 
equations 0, which can be derived by generalizing the 
Galileon theory Q and therefore is expected to be en- 
dowed with the Vainshtein mechanism. In the context of 
the Galileon, previous works focus only on the scalar-field 
equation of motion to see the profile of (the gradient of) 
(j), ignoring gravity backreaction 0, [l^.^ Since our ap- 
proach follows the cosmological perturbation theory on 
subhorizon scales the effect of metric perturbations 
can naturally be taken into account consistently. The re- 
sults in this paper can be applied to various aspects of 
cosmology and astrophysics. 

This paper is organized as follows. We define the the- 
ory we consider and then present the equations governing 
the background cosmological dynamics in the next sec- 
tion. In Sec. Ill we derive the perturbation equations 
with relevant nonlinear contributions using the subhori- 
zon approximation. We then explore spherically symmet- 
ric solutions of the perturbation equations in Sec. IV. In 
Sec. V we present some simple applications of our results 
and finally we conclude in Sec. VI. In Appendix A, we 
summarize the definitions of coefficients in the equations 
in the main text. In Appendix B, we discuss a possible 
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variety of solutions of the key equation ((50|) . In Appendix 
C, we present the Fourier transform of the perturbation 
equations. 



II. COSMOLOGY IN THE MOST GENERAL 
SCALAR-TENSOR THEORY 



We consider a theory whose action is given by 



= y d'^x^/^iCoG + An) , 



(1) 



where 



-Cgg 



Ki(ly,X)-G3{^,X)a^ 

+G4(0, X)R + Gix - {V^V,cj)f] 

-3n0(V^V,0)2 + 2(V^V,,^)3] ,(2) 



with four arbitrary functions, i^, G3,G4, and G5, of 
and X := — (9(/))^/2. Here Gix stands for dOi/dX, and 
hereafter we will use such a notation without stating so. 
The Lagrangian Cqq is a mixture of the gravitational 
and scalar-field portions, as the Ricci scalar R and the 
Einstein tensor G^i^ are included. Note in particular that 
if a constant piece is present in G4 then it gives rise to the 
Einstein-Hilbert term. We assume that matter, described 
by Cm, is minimally coupled to gravity. 

The Lagrangian ([2]) gives the most general scalar- 
tensor theory with second-order field equations in four 
dimensions. The most general theory was constructed for 
the first time by Horndeski Q in a different form than ^ , 
and later it was rediscovered by Deffayet et al. @ as a 
generalization of the Galileon. The equivalence of the two 
expressions is shown by the authors of Ref . • In this 
paper, we employ the Galilcon-like expression ^ since 
it is probably more useful than its original form when 
discussing the Vainshtcin mechanism. The gravitational 
and scalar-field equations can be found in the Appendix 
of Ref. [3. 

We now replicate the cosmological background equa- 
tions in the theory (P) [H [T!]. For = and the 
background metric ds^ = —dt^ + a^(<)dx^, the gravita- 
tional field equations are 



£ 
V 



0, 



(3) 
(4) 



where 
£ := 



V 



2XKx -K + 6X(j>HG3x - 2XG30 
-6H^G4 + 2AH'^X{Gix + XG^xx) 
^UHX^Gi^x - &H4>Gi^ 
+2i/3X(/,(5G5x + 2XG5xx) 
—6H^X (3G50 + 2XGz^x) , 

K - 2X (g30 + '(i)Gj,x) +2(iH^ + 2H]G. 
— \^H^XG4^x — 4:11 XG4X 



(5) 



(0 + 2i?0) G. 



40 



+4XG400 + AX 1^4,- 2H<pj Gi^x 
-2X (2H^(P + 2HH(P + SH^if) G^x 



-AH'X'cl)G5xx + ^HX [X - HXj G^^x 

+2 [2 iHXy+ 3H^X] G50 + 4i7X0G500, (6) 

and pm is the (nonrelativistic) matter energy density, 
while the scalar-field equation of motion is 



(7) 



where 
J : 



^Kx + 6HXG3X - 20G30 
+6H^ {Gix + 'iXGixx) - l2HXGi^x 
+2i?3x (3G5X + 2XG5XX) 
-6i720(G50 + XG50x), (8) 
— 2X ^G300 + ijiG^^x 

+6 (2i72 + ij^ G40 + QH[x + 2HX^ Gi^x 

-6H^XG5^^ + 2H^X^G5^x ■ (9) 



An overdot denotes differentiation with respect to t and 
H = d/a. 



III. PERTURBATIONS WITH RELEVANT 
NONLINEARITIES 

We work in the Newtonian gauge, in which the per- 
turbed metric is written as 

ds2 = -(l + 2$)d^2 + a2(l-2^')dx^ (10) 

with the perturbed scalar field and matter energy density. 



Pn 



0(i) + ,50(t,x), 
p„,(t)[l + '5(t,x)]. 



It will be convenient to use 

Q ■.= H- 



(11) 
(12) 



(13) 
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which is dimensionlcss. 

We wish to know the behavior of the gravitational and 
scalar fields on subhorizon scales sourced by a nonrel- 
ativistic matter overdensity S. To do so, we may ig- 
nore time derivatives in the field equations, while keep- 
ing spatial derivatives. We assume that $, 5', and Q 
are small, but nevertheless we do not neglect terms that 
are schematically written as (d^e)^ and {d^e)^, where d 
represents a spatial derivative and e is any of ^t, and 
Q. This is because L'^{t)d^e could be larger than 0(1) 
below certain scales, where L{t) is a typical length scale 
associated with G,; which may be as large as the present 
Hubble radius. [As we will see, 0{{d^e)'^) terms do not 
appear.] 

The traceless part of the gravitational field equations 
is given by 



-Y- 



where we defined a derivative operator 
1 



(14) 



(15) 



and 



-Q^^^V/Q + 2V^Qd,dkQd^d''Q 
-2dAQd^diQd''d'Q 



(16) 
(17) 

(18) 



The coefficients such as J-t-, Ai, _Bi, ... that appear in 
the field equations here and hereafter are defined in Ap- 
pendix A. 

Applying the operator djd^ to the above quantities, we 
find 



6 
1 



(19) 



(V^^V^Q - d.dj^d'd^Q) , (20) 

(21) 



where Q^^) (y^QY - (d^djQf . Thus, from the trace- 
less equation ([14| we obtain 



(J-T* - ^T* - AiQ) = 



Bi 



2a^H- 



■Q 



(2) 



^3 (v^^V^Q - d^dj<i>d'd^Q) . (22) 



The (00) component of the gravitational field equations 
reads 



Q(2) 



Ci 



(v2*v2Q-a,a,*a,9jQ) 



Q 



(3) 



(23) 



3a4i74 
where 

Q(3) := {y^qf _ 3v2q {d,d,Qf + 2 {d.djQf . (24) 
Finally, the equation of motion for (j) reduces to 

Bo 



Bi 



Q 



(2) 



Co ^(3) 



^ 0, 



(25) 



Q(2)y2^ _ 2\/^Qd,djQd'd^^ 
+2didjQd^d''Qdkd'<^. (26) 



where 



Equations ([22|) . (|23|) . and ([25|) . supplemented with the 
matter equations of motion (see Sec. IV Ap govern the 
(quasi)static behavior of the gravitational potentials and 
the scalar field on subhorizon scales. 

Note that in deriving Eqs. (|23p and we have ne- 
glected the "mass terms" {d£ /d(l))S(j) and {dS / d(j))5(f). 
These contributions could be larger than the higher spa- 
tial derivative terms, and in that case the fluctuation Scf) 
will not be excited. We do not consider this rather trivial 
situation and focus on the case where d£ /dcf) and dS/dcf) 
can safely be ignored. Though restricted to the linear 
analysis, these terms have been considered in Ref. [l3 |. 

Let us end this section with a short remark. One may 
notice that all the terms (except 5) in Eqs. ((22)) . (|23p . 
and (j25p can be written as total divergences, as 

Z^(3) = diid'^Q^^'^ - 2dj<^y^Qd'd'Q 

+2d^^dkdjQd^d'Qy 

Therefore, those equations can be integrated over a spa- 
tial domain V, and then one is left with the boundary 
terms and the enclosed mass. 



5M = p„i(<) / 6{t,^')d?x', 



(27) 



as a consequence of neglecting the "mass terms" men- 
tioned above. This fact will be used explicitly in the 
next section. 
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IV. SPHERICALLY SYMMETRIC The spherical symmetry ahows us to write 

CONFIGURATIONS 



We now want to consider a sphericahy symmetric over- 
density on a cosmological background. For this purpose 
it is convenient to use the coordinate r = ait)^ 5ijX'^x^ . 
We are primarily interested in scales much smaller than 
the horizon radius, rH <^ 1. Under this circumstance 
the background metric may be written as ds^ ~ — d<^ + 
dr^ _l_ r^dO^, where dfi^ is the line element of the unit 
two-sphere. 



where a prime denotes differentiation with respect to r. 
The gravitational field equations and scalar-field equa- 
tion of motion can then be integrated once, leading to 



2 ^' $' g' 

Ch "1 — 

r r r 

^ Q' 

\- a2 — 

r r 

Q' ^' 

ao ai a2 — 

r r r 



m\r ) m r r ' 
1 5M{t,r) _P2_fQ^ 

2 



2 



/33 *'Q' 2 71 /Q' 



3H^ 



ff2 



i72 



^2 ^ ^4 



r \ r 



(28) 
(29) 

(30) 



where we defined 



G4X ~ 0, G5 



6M{t,r) =^ 4TTp,^{t) J 6{t,r)r' dr', 
J^t/Gt, and dimensionless coefficients 



A(i) := ^, 



7.(t) :=^. 

b'T 



(31) 



(32) 



Note that Ch is the propagation speed of gravitational 
waves which may be different from 1 in general [l^ . Note 
also that in deriving Eqs. p8p -(|30 l) we have set the in- 
tegration constants to be zero, requiring that $' = vj/' = 
Q' = is a solution if SAI = 0. 

For sufficiently large r, we may neglect all the nonlinear 
terms in the above equations. The solution to the linear 
equations is given by 



clao 



ao + (2ai + c1a2)a2 ' 

ao + aia2 ^ 
ao + (2ai -I- c\a2)a2 ' 



ai -I- Cf^a2 



ao + (2ai + c\a2)a2 



(33) 
(34) 

(35) 



where wc defined /i SM/SttGt- In this regime, the 
parametrized post-Newtonian parameter 7 is given by 



ao + aia2 



c^ao 



a 



2 ' 



(36) 



which in general differs from unity. 



A simple example for which nonlinear terms can op- 
erate is the model with G4X = = and G^x 0> 
i.e., 



C = Gi{<t>)R + K{<t>,X) - G:i{4>,X)Oq 



In this case, we have Pi = h = - 
c\ = l. We also have the relation 

/3o = Y + «2(^ 0). 



70 



(37) 

71 = and 



(38) 



The Lagrangian (|37[) corresponds to a nonminimally 
coupled version of kinetic gravity braiding [l5| , and has 
been studied extensively in the context of inflation [l6l | 
and dark energy/modified gravity [T7l-[l9j. Of course the 
nonminimal coupling can be undone by performing a con- 
formal transformation, but in the present analysis we 
have no particular reason to do so. Note that even in 
the case of G4 =const the scalar is coupled to the cur- 
vature at the level of the field equation, which signals 
"braiding." 

Using Eqs. (EH) and (HH, Eq. ^ reduces to a 
quadratic equation 



where 
B := 



B_ 

4;9o 



Q' 



2Q' 



ao H- 2aia2 



9 ' 



C 



ai + a2 



ao -|- 2aia2 



(39) 



(40) 
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Equation ((39)) can easily be solved to give 



At short distances, <^ BC^/H'^, one finds 



(41) 



(42) 



In order for this solution to be real, we require 

BC>0 ^ Gsx {XG-sx + G40) > 0. (43) 
In this regime, the metric potentials are given by 



GnSM H{ai+a2) 2BCGnSM ^^^^ 
$ ^ -\l , (44) 



B 



GnSM Ha2 I2BCGnSM 



where 



264 



(45) 



(46) 



If BC - 0(1), the typical length scale n {fi/H'^y/^ 
can be estimated using Eq. p6|) as 



1/3 



120 



H 



HoV^' fSM\ 



1/3 



pc. 



(47) 



where Hq = 70km/s/Mpc. 

Note that in this case the Friedmann equation can be 
written as 



3iJ^ = SttGcos (Pm + P4>) ' 



(48) 



where Gcos = l/167rG4 = Gn and P4, = 2XKx — K + 
6X4>HG3x -2XG34,-6H<j)G4^. Thus, the effective grav- 
itational coupling governing short-distance gravity is the 
same as the one in the Friedmann equation. 

The Vainshtein mechanism successfully screens the ef- 
fect of the fluctuation Scj), so that the two metric po- 
tentials coincide and exhibit the Newtonian behavior at 
leading order. However, Gat is time-dependent since it is 
a function of the time-dependent field 0(t), which means 
that the Vainshtein mechanism cannot suppress the time 
variation of Gjv in a cosmological background. This fact 
was first noticed in Ref. [13] ■ We will discuss this point 
further in the next subsection. 



B. 



G5X = 



With G5X = 0, the problem reduces to solving a cu- 
bic equation one can handle. Indeed, using Eqs. ([25)1 
and (1211) to remove <&' and from Eq. ([30]), we obtain 
the cubic equation for Q': 



H'^CaP 



0, 



where r-independent coefficients are defined as 



(50) 



c — 
Cp := 
Ci := 



213,132 +cil3i- jo' 



f3i + cl(32 



2/3,(32 + cf^l3l - V 
gp + (2qi + c\a2)a2 
2/3i/32 + c2/3| - 70 ■ 



^ _ 2/?o + 3 (ai/32 + a2/3i + c^a2/?2) 

2 (2/31/32 + 4/31-70) ■ ^ ' 

We note the expressions for ^P' and $' in terms of Q': 

^' = - (ai + 4"2)g' - (/3i + 4/32) ^" 



*' = 4-"2Q'-/32 



(52) 



Linearizing Eq. (|5(I)) at ^ {Cp/Ci)p/ H"^ , one obtains 
the solution ([551) as expected. This will be matched to 
one of the following three solutions at short distances: 



Q' 



-hJc 



c 



(53) 



If simply Ca ^ Ci ^ C2 ~ 0{1),^ the two regimes 

arc connected at around r ~ r* , where 



Cl P 



1/3 



for Q'~±HJCp 



P 



{-TrTHj for g'^- — — . (54) 



Ci 



If C/3 > and Ci > 0, the solution with the bound- 
ary condition (j35p at large r can be matched cither to 
Q' ~ +H^Cfjp/r or to Q' ~ -Hy/Cpp/r. We caU 
this situation Case I. This is possible for {C2,Ca) in the 
shaded region in Fig. [TJ ForC^ > (respectively Cq, < 0), 
the short-distance solution is given by +H ^J'CjjJxJr (re- 
spectively ~H ^JCppJr). Outside this region one cannot 



Let us consider a class of models with G^x = 0. In 
this case, we see that /33 = 71 = 0. For the other nonzero 
coefficients we have the following relations: 

4 = l + 2/3i(^l), /3i -f /32 + 270 = 0. (49) 



^ This is probably the most natural case if one considers a model 
that accounts for the present cosmic acceleration and the coeffi- 
cients are evaluated at present time, because in that case there 
is only one typical length scale L = Hq^ . 
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FIG. 1: Relation between the coefficients and the short- 
distance solution for Case I. In the shaded region one gets 
a real solution whose behavior at short distances is noted. 
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7 # 1 (Eq. 36) 
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FIG. 3: Coefficients for which a real solution exists [Case IF). 





C2 = -V^ 












c„ n'-'r 
















Cg < 0, Ci > 



FIG. 4: Coefficients for which a real solution exists {Case IT). 



FIG. 2: r^$', r^*', and r^Q' as a function of r. H'^ = 1. where 
cl = 1.2, Qfi = 02 = 1, Qo = 2, ^32 = 2, /X = 10"^°. Inside the 
Vainshtein radius 7 = 1 is reproduced, but the concrete value 
of the radius depends on the model under consideration (see 
Eqs. (HZll and (pl)). 



find a solution that is real for r G (0,oo). A typical 
behavior of the Case I solution is plotted in Fig. [2] If 
CpCi < 0, then the solution can be matched only to 
Q' ~ -{Ca/Cp)H^r/2. We call this situation Case II. 
This is possible for (C2,Cq) in the shaded region in Figs. [3] 
{Cfs > 0, Cl < 0) and|l(C/3 < 0, Ci > 0). Outside this 
region no real solutions can be found. If C/3 < and 
Cl < then no real solutions can be found either. 

Let us then evaluate the metric perturbations for each 
solution Q'. We begin with the Case I, Q' ~ ±y/C^JI/r. 
The metric potentials at short distances are given by 



5M 



SM 
8t:Qt t'^ ' 



(55) 



2/3i/32+c2/3|-7o' 

/^ife - 70 
2/?i/32 + c2/3|-7o- 



(56) 
(57) 



Although the coefficients look apparently different, now 
we use the relations (l49l) for the first time to show that 



(58) 



0$ — Oxp — 



i.e., the two metric potentials actually coincide. We thus 
obtain the Newtonian behavior 



SttQt 



(> 0). (59) 



It is interesting to note that the above conclusion holds 
even for the generic propagation speed of gravitational 
waves, c| 7^ 1. Explicitly, one finds 



1 



N 



2 (G4 - 4XGix - 4A2G4XX + SXG^^) 



(60) 
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As in the case of the previous subsection, Gn is in gen- 
eral time-dependent, as it is a function of time-dependent 
and X — (j)'^{t)/2. We thus illustrate how the Vain- 
shtein mechanism fails to suppress the time variation of 
Gjv in a cosmological background within the context of 
some generic scalar-tensor theories minimally coupled to 
matter. The claim was originally suggested using the 
Einstein frame action in Ref. [20|. Here we explicitly 
give the concrete formula with which one can evaluate 
the time variation of Gn for a given model. 

If Gat happens to vary very slowly, we can say that 
the usual Newtonian gravity is reproduced in the vicin- 
ity of the source. However, in general, one expects that 
Gn varies on cosmological time scales. The time varia- 
tion |G/G| is constrained from lunar laser ranging exper- 
iments to be \Gn/Gn\ < 0.02i?o HH. 

At this stage it is interesting to look at the background 
evolution for the models with G^x = 0. The (modified) 
Friedmann equation ([3]) can be written as 



For the solution Q' ~ -{Ca/Cp)H'^r /2, we find 



(61) 



where the gravitational coupling in the Friedmann equa- 
tion read off from the above exactly coincides with the 
expression for Gn, 



G(-ns — G 



N, 



(62) 



and 



2XKx -K- 2XG30 
+6H (x^G3x - 2X0G40X 



The situation here is the same as what we have seen 
in the previous subsection. We refer to a constraint in 
Ref. i22j, obtained by translating the big bang nucle- 
osynthesis (BBN) bound on extra relativistic degrees of 
freedom, as 



1 



G 



W BBN 



G 



I now 



^0.1, (63) 
is evaluated at the 



where G^vIbbn (respectively, G^vlnow 
time of BBN (respectively, today). 

Having thus seen that the Newtonian behavior is re- 
produced with time-dependent Gn, let us then evaluate 
leading order corrections to the potentials. In this case 
we need to keep the subleading term in Q': 



Q' ~ ±HJCp^ 



4Cs 



From this we obtain the corrections A$' 
GnSM/v^ and A^-' = - GnSM/v^ as 



(64) 

= $' - 



A*' = tH 



"2 + 7^ (Cq - 2C/3C2) 



2 , /3l + 4^2 

ai + c^a2 H ;r;; (Ca 

ZLr 



2GnSM 



2GnSM 



(65) 



(66) 



SM 



+ 0(r), ^' 



1 m 



0(r), (67) 



implying that the parametrized post-Newtonian param- 
eter 7 is given by 7 = l/c|. Therefore, cf^ is tightly con- 
strained from solar-system tests in this case: |1 — 7I < 

2.3 X 10-5 m. 

When the coefficients Cq, Cp, Ci, and C2 have hi- 
erarchies in their values, we find a variety of solutions 
to Eq. ([50]) on an intermediate scale between the linear 
regime at large r and the small r limit of Case I or Case 
II. The details are summarized in Appendix B, which 
could be potentially confronted with observations. 



C. Gsx / 

Let us finally discuss the most general case where all 
the coefficients in Eqs. are nonzero. Although 

one can still eliminate $' and ^' to get an equation solely 
in terms of Q' , the resulting equation will be a sextic 
equation. This hinders us from analyzing a variety of 
possible solutions in detail. However, for /Sa 7^ and 
7i ^ one can show that there is no solution such that 
~ V]/' ^ l/r^ on sufficiently small scales. To show 
this, one substitutes $' ~ ^f' l/r^ to Eq. dH]). The 
second term in the right-hand side can be compensated 
by the other provided that Q' ^ r or Q' ^ 1/r^. If Q' ^ 
r, one cannot find a term that compensates the forth 
term in the right-hand side of Eq. (PH)) which behaves as 
^'vj/'/r^ ^ 1/7.6 Qi ^ then one cannot find a 

term that compensates the last term in the right-hand 
side of Eq. (|29l) . Thus, there is no consistent solution 
with ~ V]/' 1/7-2 Qj-^ sufficiently small scales. This 
implies that the typical length scale associated with B3 
and Gi must be as small as 0(100 ^m) [IJ], though it 
is uncertain whether or not we can have the Newtonian 
behavior of the potentials on intermediate scales. 



V. APPLICATIONS 

A. Evolution of density perturbations 

Since it is assumed that matter is minimally coupled 
to (f), no modification is made for the energy conservation 
equation and the Euler equation for matter. Therefore, 
the nonlinear evolution equation for S is given by 



• 4 <j2 
S + 2HS--- 



3 1 



V2 

= (1+5) $. 



(68) 



However, as we have seen in the previous section, the 
relation between $ and S is modified. One may tackle 
the nonlinear equations using the perturbative approach 
(e.g., [IH). We derive the Fourier transform of the non- 
linear equations in Appendix C. 
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0.001 0.01 0.1 1 10 100 



r (/i^'Mpc) 



FIG. 5: The ratio of the circular speed V'^ of a test particle 
in the generalized model to that in Newtonian gravity, as a 
function of radius r(ft~^Mpc). We used the Hubble parameter 
at present H = 100 h km/s/Mpc where h — 0.7, and the 
velocity dispersion a = 100 km/s. Each line corresponds to 
the following cases: (a) c| = 1, qq — —2.0, 02 = 1.0, /So — 1.0, 
ai = /?i = /?2 = 70 = 0.0, (6) cl = 0.9, ao = -2.0, ai = 0.1, 
02 = 0.9, Po = 0.9, /?i = -0.05, /32 = 0.1, 70 = -0.025, (c) 
c| = 0.9, ao = 2.0, ai = 1.0, 0^2 = 0.1, fio = 1.0, /3i = -0.05, 
132 = 0.1, 70 = —0.025. Dotted and dashed lines are given by 
Eq. ([73) and (O, respectively. 



For simplicity, we here consider the cases G4X = = 
G5 and G5X = as described in llV Al and llVBI Although 
the mass is proportional to r, one can check that three 
solutions at short distances remain the same as Eq. (|53p . 
To see the effects of modification of gravity, we consider 
the velocity V{r) of a test particle in a circular motion 
with radius r, V'^{r) — r$', which reduces to 2a^ in New- 
tonian gravity. In the generalized model, V'^{r) depends 
on the radius r, whose asymptotic behavior can be found, 

V - (73) 

for large r, and 

~ (74) 

for small r, where ~ 2a^. A typical behavior of 
the circular speed divided by the one in general rela- 
tivity, V^/V^, is demonstrated in Fig. [S] The line (a) 
represents the minimally coupled model, corresponds to 
G4 = MlJ2 and G5 = 0. The lines (6) and (c) show 
the model with G^x = and the parameters in (h) and 
(c) are chosen so that the solutions at short distances 
become Q' ~ -{Ca/Cp)H'^r /2 and Q' ~ +H^Cf)fj./r, 
respectively. As one can see in Fig. [Sj the Vainshtein 
radius ^ 0{1) Mpc for cr = lOOkm/s. 



For spherical perturbations, V^$/a^ can be expressed 
in terms of S using the results in the previous section. It 
follows that 



(69) 



where the effective gravitational coupling for large-scale 
perturbations (but well inside the Hubble horizon) is 



G, 



1 



eff • = 



8nQT ao + {2ai + c^a2)a2 ' 



and 



V2 

— $ -> AnGNPmS 



(70) 



(71) 



for small-scale ones. 



B. Halo Model 

We consider a simple halo model to investigate a char- 
acteristic feature of general second-order scalar-tensor 
theories [26j . For simplicity, let us assume the density 
of matter follows 



Sp{r) 



27rGjvr2 ' 



(72) 



where a is the parameter of the velocity dispersion. In 
this model we have SM{r) — 2a^r/GN- This model is 
the singular isothermal sphere in Newtonian gravity. 



VI. CONCLUSION 

Based on the most general scalar-tensor theory with 
second-order field equations, we have studied metric per- 
turbations on a cosmological background under the in- 
fluence of the Vainshtein screening mechanism. We have 
derived the perturbation equations with relevant nonlin- 
earities by taking into account the effects of cosmological 
background. We have clarified how the Vainshtein mech- 
anism operates in the two subclasses: (i) G^x = G5 = 
and (ii) Gzx ~ 0. The situation in the first case 
Gix = G5 = is very similar to the Vainshtein mech- 
anism in the Galileon theory, which contains only XDcj) 
in the Lagrangian. However, the second case G^x = is 
considered for the first time in a cosmological background 
in this paper. We have explicitly shown that below 
the Vainshtein scale there are three possible solutions 
for Q': Q' ~ ±Hy/t^^ and Q' ~ -{Ca/Cf})H'^r /2. 
We have found that two metric perturbations coincide 
well inside the Vainshtein radius in the first case, 
while the parametrized post-Newtonian parameter 7 is 
related to the propagation speed of gravitational waves, 
c^, well inside in the second case. In both cases, New- 
ton's constant Gn, its time variation \Gn /Gn\j and the 
parametrized post-Newtonian parameter 7 can be con- 
strained from BBN and the experiments such as lunar 
laser ranging. These could provide powerful constraints 
on the most general scalar-tensor theories. In the case 
G^x 7^ 0, we have demonstrated that the inverse-square 
law cannot be recovered at sufficiently small r. 
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Appendix A: Coefficients in the field equations 

Here we summarize the definitions of the coefficients in the field equations: 

Gt 

e 

A2 

Bo 



2 6*4 — 2XGix — X (^H(f)Gzx — G^^^ 

-4>XG3x + 2iJG4 - 8HXG4X - SHX^Gixx + (I^Gi^ + 2X(PG44,x 
{5XG5X + 2X^G5xx) + 2HX (SGs^ + 2XG50x) , 

e e Jt £+r 



1 dGr 
H dt 

X 



H 



+ 3 (X + 2HX) Gixx + 2XXGixxx - 'i(l>Gi^x + 2(j)XGi^xx 



2HX 



H + Hn X 



H + H^j <j>G5x H 
— (/)XG500x — X (^X — 2HXj G^cjyxx 
2X G4X + (f) {G5X + XG^xx) — G,50 + XG^^x 
—2X (^G^x + 2XGixx + H4)G5x + H4)XGzxx — G50 — XG^^x 
H^XG^x, 



G5XX + H<t>XXG^xxx -2[X + 2HXj 



G: 



54>X 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 
(A6) 



2X^G^ 



2X^ 



AXX 



2>i>G^xx + 'i>XG^xxx — 2G^(j>x + XG^^xx 



H^XiGsx+XG^xx). 



(A7) 

(A8) 
(A9) 
(AlO) 
(All) 
(A12) 



Appendix B: Intermediate regime for Gsx = 

Here we would like to point out that there could be an interesting intermediate regime where the quadratic term 
in Eq. ([SH]) comes into play so that we have the solution 



This regime can be found in the range 



C£JJ_ 



1/3 



1/3 



CaC2 M 



1/3 



1/3 



(Bl) 



(Case I), 



(Case II). 



(B2) 
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This intermediate regime can be seen if Ca or C2 {Cp or Ci) is sufficiently large (small) compared with others. It is 
also interesting to see the behavior of the metric perturbations in the intermediate regime (|B2[) . In this regime, the 
metric potentials can be obtained by substituting Eq. (jBl| into the relations (|52|) . 



ii 



(B3) 
(B4) 



The metric potentials in this regime differ from Eqs. (|59p and (j67p . and the parametrized post- Newtonian parameter 



7 



2C2Cj3 — P2CaC^ 



2clC2Cp ~ CI/32C2C0. - P1C2C0, ' 



(B5) 



is not equal to unity. 

We also notice another intermediate regime if Ca is sufficiently large compared with the other coefficients. In this 
regime Q' becomes constant, 

1/3 



Q' 



This solution can be seen between Eq. (jBip and Eq. ((53 

1/3 



\CI H^j 

The metric potentials in this regime are given by 



Cg ^1 \ 



1/3 



•^h ^2 ■ 



(B6) 



(B7) 



(B8) 



Thus, the parametrized post-Newtonian parameter is 7 = in this regime. If these regimes (|B2p and (|B7[) include 
our solar-system scales, it is possible to constrain the parametrized post-Newtonian parameter 7 as in the former case. 



Appendix C: Equations in the Fourier space 

In this Appendix we summarize the coupled equations for the evolution of the matter density perturbations in 
Fourier space, which will be useful in the perturbative approach (25| . The matter density perturbations follow 

|J + -V-[(l + <5)v]=0, (CI) 
ot a 

_Z + i/v + -(vV)v = V$, (C2) 

ot a a 

where v is the velocity field. Assuming the irrotational fluid, we introduce the velocity divergence = V • v/(ai7). 
Then, due to the Fourier transform, the above equations can be rephrased as 



H dt 



P) = -j^ J dk,dk2S^^\k, + k2 - p) (^1 + ^jP-^ e(k2)<5(ki), (C3) 



1 d9(p) I h\ ^, , ^, . 

+ 2 + T7J^(p)-7lW*(p) 



H dt \ / a2^2 

= / rfk,dk,^(3)(k, + k2 - p) ( ('^^ - '^^^jj^^^ + '^^l' ) eik,)e[k2). (C4) 

Similarly, Eqs. a.nA ^ give 

' (J-T*(p) - Qr-^iv) - ^iQ(p)) - 2a^g2 (2^)3 / dk,dk25^^\k, + kz - p) {klkl - (ki • ^2?) Q{k,)Q{k2) 

^ J dkirfk2,5(3)(ki + k2 - p) {kjk^ - (ki • k2)2) Q(ki)$(k2), (C5) 



a2iJ2 (2 
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B9 



2a^H^ (2 
S3 



y" dkidk2^(3)(kl + k2 - p) (fc2fc2 „ (ki . k2)2) Q(ki)Q(k2) 

a2iJ2 (2"^ / rfkidk2<5(3)(ki+k2-p) (fc2A:|-(ki.k2)2)Q(ki)vI/(k2) 

^ y" dkidk2dk3(5(3) (ki + k2 + kg - p) 



xo(ki)g(k2)g(k3) 



-kfklk^ + 3fc?(k2 • k3)2 - 2(ki • k2)(k2 • k3)(k3 • ki) 



(C6) 



1 



(27r)3 
B2 1 
a2iJ2 (27r)3 

^3 1 
a2ij2 (27r)3 

Co 



So 1 


J dkidk2(5(^Hki +k2 - 


2i72 (27r)3 


P) (fc?fc| - 


(ki •k2)2)«'(ki)Q(k2) 


P) {klkl - 


(ki •k2)2)$(ki)Q(k2) 


P) - 


(ki •k2)2)$(ki)4'(k2) 



J rfkidk2rfk3<S(3)(ki+k2+k3-p) 

klk^kl + 3fc?(k2 • kaf - 2(ki • k2)(k2 • k3)(k3 • ki) 

^ J dkidk2rfk3(5(3nki+k2+k3-p) 

-klk^kl + (ki • k2)^kl + 2fc2(k2 • kg)' - 2(ki • k2)(k2 • k3)(k3 • ki) 



a^H^ (2 

xQ(ki)Q(k2)g(k3) 

a4iJ4 (2 

xQ(ki)Q(k2)$(k3) 



= 0, 



(C7) 



respectively. 
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